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ABSTRACT 

We compare the measured peculiar velocities of 98 local (< 150 h^ 1 Mpc) type la su- 
pernovae with predictions derived from the PSCz survey. There is excellent agreement 
between the two datasets with a best fit /3j (= fl^f/bi) of 0.55 ± 0.06. Subsets of the 
supernovae dataset are further analysed and the above result is found to be robust 
with respect to culls by distance, host-galaxy extinction and to the choice of reference 
frame in which the analysis is carried out. Alternative methods of determining j3i in- 
cluding density-density comparisons, dipole measurements and WMAP-based results 
are also discussed. We conclude that most recent determinations are consistent with 
a value of /3j = 0.5. 

Key words: cosmology: observations - cosmological parameters - dark matter - 
galaxies: distances and redshifts - large-scale structure of Universe 



1 INTRODUCTION 

Peculiar motion studies are a powerful tool for examining 
the underlying mass distribution of the local universe. In 
the linear regime, where density fluctuations are small, the 
mass over-density, S m , can be related to the fluctuation in 
galaxy number- density, S g , via S g (r) = bS m (r) where b is the 
linear bias parameter. This bias parameter together with the 
cosmological mass density parameter, fl m , can be used to 
predict peculiar velocity fields from all-sky redshift surveys 
via the equation 



19991) . ISchmoldt et alJ il99Sh and iRowan-Robinson et alJ 
2000). Recent comparisons of these fields with peculiar ve- 
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The dimensionless quantity f3 (= f2^ 6 /6) scales linearly with 
the pre dicted velociti es and is the only free parameter of the 
model JPeebleslll980l) . Hence /3 can be directly determined 
from the comparison of measured peculiar motions with pre- 
dictions made from galaxy density fields. 

All-sky galaxy samples derived from IRAS satellite data 
have been used extensively to map the local density field. 
Currently the most compl ete redshift survey of IR AS sources 
is provided by the PSCz jSaunders et al]l2000bh . This sur- 
vey consists of redshifts for 15,411 galaxies uniformly dis- 
tributed over 84.1% of the sky with a median redshift of 
8500 km b -1 . The PSCz survey's depth, excellent sky cov- 
erage and density allow for the reliable mapping of the dis- 
tribution of galaxies in the local universe. Several indepen- 
dent determinations of the PSCz density a nd velocity fields 
have therefore been made; most notably bv lBranchini et ail 



locity measu rements typica lly yield values of f3i in the range 
0.4 - 0.6 fsee lZaroubill2002l) . 

A significant source of error in determining (3 arises 
from the uncertainty in the peculiar velocity measurements. 
Galaxy distance estimates from the Tully-Fisher and Fun- 
damental Plane relations are subject to errors that are 
typically ~ 20% per galaxy. At depths greater than ~ 
50 h~ l Mpc this is considerably larger than the peculiar 
velocities of the individual galaxies. With distance errors 
less than 10%, Type la supernovae (SN Ia) are less susc epti- 
ble to inhomogeneous Malmquist bias jHudsor]ll994ah and 
hence offer an important alternative probe of the local ve- 
loc ity field. An early attempt to use SNIa was carried out 
bv lRiess et alJ il997t) who compared the peculiar velocities 
of 24 SNIa with the vel ocity fields predict ed from the 1.2 Jy 
IRAS redshift su rvey l|F isher et al.lll995|) and the Optica l 
Redshift Survey l|Santiago et al] 1199.4 iBaker et alJll998n . 
They derived /3j = 0.4 ± 0.15 and j3 = 0.3 ± 0.1 respec- 
tively, with the relatively large error resulting from the small 
sample size. 

iTonrv et all (12003 



have recently produced a homoge- 
nized compendium of 230 SNIa for constraining cosmological 
quantities. The release of this compendium presents a new 
opportunity to measure f3 with a significantly smaller error. 
In this paper we c ompare pecu l iar ve locities measured by the 
local SNIa in the lTonrv et alJ i2003T) sample to the peculiar 
velocity field deri ved from the smoothed PSCz density field 
as determined bv lBranchini et all Jl99Sl) . 
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We briefly describe the lTonrv et al] J2003f) compilation 
of SNIa in Section [5] and review the derivation of the PSCz 
velocity field in Section |H] In Section 2] we describe the 
derivation of fti. The robustness of this result is analysed 
in Section and finally, in Section |S| we discuss the results 
and present our conclusions. 



2 THE SNIA DATASET 

The iTonrv et al] (|2003) dataset is a homogenized com- 
pendium of 230 SNIa compiled f r om many recent s tudies . 
Most notably from t helJhal ]2002l) Iperlmutter et al.Ul999t). 
iHamuv et al j (ll996l ). lRiess et al](ll999r ) and lGermanv et alJ 
(2004) datasets, which comprise the majority of the data. 
Using a variety of fitting techniques such as MLCS (Riess et 
al. 1998 and the work of Jha and c ollaborators) and dml5 
jGermanvll200ll) . lTonrv et al] (120031) have re-calculated the 
relative SNIa distances where the original photometric data 
is available. The systematic offsets of each dataset were re- 
duced by minimising the differences between all pairs of 
datasets where overlaps exist. The residuals of this fitting 
procedure are 0.02 mag or better for th e majority of the 
samples. Table 15 of ITonrv eta!] (120031) lists the redshift 
(logos), luminosity distance (logdHo), distance error and 
hos t galaxy l/-band ext inction (Ay) for each SNIa. 

ITonrv et all (I2003T) fix the zero point of nearby SNIa 
(0.01 < z < 0.1) by assuming an 'empty universe' (Q m = 
0, f^ A = 0) cosmology. For our analysis, we have converted 
the ITonrv et al] (120031) quoted distances to a ACDM cos- 
mology (57a = 0.7, £l m = 0.3). However, the derived fti is 
unaffected by the choice of cosmology. 

In this paper we only consider the 107 SNIa that 
lie within 150 hT 1 Mpc as the PSCz density field is 
incomplete at greater distances for all galactic latitudes 
jBranchini et al.lll999jl . We further restrict the sample to 
SNIa with extinctions Ay < 1.0 mags, for reasons discussed 
below. These selection criteria leave 98 SNIa, which we refer 
to as the "default sample". The median distance error for 
this local SNIa sample is ~ 8%. 



3 THE PSCZ VELOCITY FIELD MODEL 

iBranchini et al] (I1999T) used the PSCz redshift survey to 
determine the density and peculiar velocity fields in real 
space in a self-consistent way by using equation (0 un- 
der the assumption that mass follows the number den- 
sity of IRAS galaxies. These fields are smoothed with 
a Gaussian filter of radius 5 h~ Mpc. Analysis by 
iBerlind. Narayanan fc Weinberg] (l2000f) indicates that this 
smoothing radius shoul d yield unbiased result s for fti. In 
an independent analysis, ISchmoldt et al] dl999F) derived the 
PSCz velocity and density fields by using a Fourier-Bessel 
approac h. They found the resu lting fields to be consistent 
with the IBranchini et al] Jl99Sf) fields used here. 

The integral in equation Q extends over all space. The 
PSCz survey, however, does not extend to infinite depth, 
nor does it have data in the Zon e of Avoidance (ZoA). For 
the ZoA, IBranchini et al. | jl999l) ha v e imp lemented a simi- 
lar approach to that of lYahil et al] al991f) by dividing the 
region (|6| < 8°) into bins of 10° latitude by 1000 km s -1 . 



These bins are then populated with enough synthetic galax- 
ies to reflect the number density of the corresponding bins 
at greater |fe|. The systematic effect on the derived value 
of fti due to this i nterpolat i on pro cedure can be estimated 
from the results of iHudsonl (Il994lj) . He compared ft values 
derived from an optically-selected density field with a larger 
ZoA (|6| 5j 12°) using different techniques to account for the 
missing structure. Only an 8% difference was observed be- 
tween the ft value derived from the interpolated density field 
and that derived from a density field in which the ZOA was 
assumed to be at average density. Since the average density 
assumption is rather extreme, this result may be taken as 
an upper limit on the systematic uncertainty. Therefore, as 
the PSCz ZoA is only two-thirds the thickness of this ZoA, 
we might expect a systematic uncertainty on our result of 
the order 5%. This is considerably smaller than our random 
errors. 

As stated previously we have truncated the PSCz ve- 
locity field at 150 h^ 1 Mpc due to increasing shot noise. 
Sources beyond this depth, however, may still contribute to 
the LG's motion. Because the statistical weight of the SNIa 
sample is dominated by nearby objects these external con- 
tributions can be modelled as a dipole term. For peculiar ve- 
locity comparisons in the LG frame this dipole term cancels 
out as the motions of the LG and SNIa are affected in the 
same way. LG-frame comparisons assume, however, that the 
LG's motion is exactly given by linear theory. In practice, 
the LG is expected to exhibit a nonlinear 'thermal' compo- 
nent to its velocity that is not well modelled by linear theory. 
An alternative to the LG-frame comparison is to omit the 
LG from the analysis entirely. This can be achieved by fit- 
ting the SNIa peculiar velocities in the CMB frame with an 
additional dipole component to allow for contributions not 
included in the PSCz density field. Ideally, analyses in both 
these frames should produce similar results. However, due 
to the larger uncertainty in the CMB analysis, we regard 
the LG result as a more reliable solution. 



4 DETERMINING fti 

There is a very good agreement between the peculiar veloc- 
ities measured by the SNIa and predicted from the PSCz. 
This is shown in Fig. 0where the scatter around the Hubble 
flow before and after the PSCz velocities for fti — 0.5 are 
removed. In the range 20 — 80 /i -1 Mpc, where the major- 
ity of SNIa lie, the removal of the predicted PSCz peculiar 
velocities reduces the rms scatter around the Hubble flow 
from 490 km s~ x to 390 km s _1 . In Fig. □ nine SNIa with 
Ay > 1.0 are plotted as open circles, three of which are 
distinct outliers. In our analysis we have chosen to exclude 
these objects because we expect that their errors are under- 
estimated. 

To determine fti in the LG frame we minimise the \ 2 
relation: 




where Vi is the peculiar velocity of the i supernova, Vi,pscz 
is the PSCz-predicted peculiar velocity which depends on fti 
from ad is the distance error and a cz incorporates both 
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Figure 1. The Hubble flow residuals for all 107 SNIa lying within 
150 h~ 1 Mpc in the LG frame. The upper panel shows the orig- 
inal uncorrected data whilst the lower shows the data with the 
predicted PSCz peculiar velocities removed. Note the reduction in 
scatter, particularly in the distance range 20-80 h~ 1 Mpc. SNIa 
with host-galaxy extinctions Ay ^ 1.0 are plotted as open circles 
whilst filled circles show the default sample used in this paper. 



an estimate of the error in redshift determination as well as 
errors in the PSCz predictions due to shot noise or non-linear 
peculiar velocity contributions. 

Various studi es have adopted different schemes for <r cz . 
iRiess et alJ ( jl997f) , adopt a value of 200 km s -1 for all the 



use values of 150 km s 
Blake slee et alJ l|l999F) also ac- 



SNIa, whilst Blakeslee et al 
and 200 km s -1 . However 
count for the extra velocity dispersion of cluster galaxies 
using two different approaches. Their 'Trial 1' method adds 
in quadrature an extra factor of a c i(r) = croA/l + (r/ro) 2 
to cr cz where a = 700 (400) km _1 and r = 2 (1) Mpc for 
galaxies in Virgo (Fornax). Their 'Trial 2' scheme uses the 
standard (j cz but resets the individual galaxy velocities for 
group members to the group-average velocities as listed in 
iTonrv et al] il997t) for 37 separate clusters. In our analysis 
we extended both these techniques to account for galaxies 
which lie near one of the X- ray selected cluster s of the NOAO 
fundamental plane survey iSmith et al] l20041. 

Table lists the derived f3i values for these different 
weightings for our default sample. The la quoted errors are 
calculated from bootstrap re-samples of the dataset. If the 
nine Ay > 1.0 SNIa had not been removed, the resulting \ 2 
would be larger by ~ 40. 

Increasing the redshift error er cz for SNIa lying close to 
nearby clusters has a sizeable effect on the x 2 but appears to 
have no significant effect on the value of 0i. Overall, little 
variation from the preferred value of /3i = 0.55 ± 0.06 is 
observed and j3i is effectively independent of the weighting 
schemes used. 



Table 1. "Redshift error" , <r cz , comparison for the default sample 
of 98 SNIa in the LG frame. The errors have been determined 
from the la deviation in the distribution of the medians of 1000 
bootstrap re-samples. 
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Figure 2. Comparison of SNIa peculiar velocities to PSCz pre- 
dicted peculiar velocities in the range h' 1 Mpc to 150 h' 1 Mpc 
with Ay < 1.0, cr 2 z = 150 2 + cr 2 and (3 = 0.55. The top panel 
shows comparisons in the LG frame, and the bottom panel shows 
the comparison in the CMB frame (without the extra dipole com- 
ponent). The size of the data point is inversely proportional to 
the total error (a = ^/tr 2 + cr 2 z ) on each SNIa. The smallest 
and largest circles correspond to values of a = 1290 km s — 1 and 
170 km s -1 respectively. The lines indicate a 1:1 ratio. 



In order to determine j3i in the CMB frame an extra 
dipole component is added as an extra free parameter in 
the minimization of equation JS). Using the default sample 
with a cz given by 'Trial 1' as y/ 150 2 + oj, the best fit has 
Pi = 0.48 ± 0.09 and V dipo i e = 206 ± 97 km s" 1 towards 
I = 290° ± 25°, b = 0° ± 18°. This extra dipole component 
is consistent with zero but is also consistent with the value 
of V dip oie = 372 ± 127 km s" 1 towards I = 273° ± 17°, 
b = 6° ± 15° as found bv lHudson et al] (1200411 for the SMAC 
sample. The calculated value of f3i agrees well with the result 
derived in the LG frame. 

The good agreement between the observed and pre- 
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dieted peculiar velocities in both the LG and CMB frames 
is shown in Fig. [5] If the peculiar velocities predicted by the 
PSCz and observed from the SNIa are in exact agreement 
for the chosen value of Pi, the SNIa would be expected to 
lie along the 1:1 line. This trend is indeed observed. The dif- 
ferences between the measured and predicted velocities are 
as expected given the errors in both distance and velocity 
measurements, i.e. the data is consistent with a reduced xt 
of ~ 1. Thus the two datasets agree exceptionally well. 



5 ROBUSTNESS 

To assess the robustness of the derived f3i we have exam- 
ined various sub-samples of the local SNIa dataset. Unless 
otherwise stated all sub-sa mples use o ur default sample in 
the LG frame with er cz = t^/150 2 + a^i determined using the 
'Trial 1' approach. Table|5|lists the best fit j3i together with 
the associated x 2 f° r each sub-sample. 

Importantly, /3j is found to be independent of the dis- 
tance range considered. Any derivation of f3 is expected to be 
strongly weighted by the very nearby SNIa where measure- 
ment errors are smallest. Hence we have tested the depen- 
dency of our calculations on SNIa at different distances by 
dividing the data into two distance ranges. The position of 
this division is chosen such that the bootstrap errors on each 
derived Pi are of similar magnitude. For a distance range of 
— 30 /i" 1 Mpc we derive a value of Pi = 0.55 ± 0.07 and for 
30 - 150 h' 1 Mpc, fa = 0.54 ± 0.10. Tablegjalso includes a 
variety of different distance ranges all of which yield similar 
values of pi (0.49 < (3 < 0.58). 

The determination of pi is also revealed to be indepen- 
dent of the cull by host-galaxy extinction with Pi varying 
by only ±0.05 for culls down to Ay < 0.3. It is found that 
the reduced Xv is ~ 1 for all culls of host-galaxy extinc- 
tion < 1.0. Overall, for all the sub-samples considered, Pi is 
found to range by only 0.10. 

Another source of bias which we do not account for in 
our analysis is inhomogeneous Malmquist bias. Generally, 
not correcting for this will lead to higher values of p. How- 
ever, this bias scales with the square of the distance error. 
Thus for the SNIa sample used here we expect this bias to 
be considerably smaller than the random error in Pi. 



6 DISCUSSION 

Tabled lists a representative set of recent determinations of 
pi from comparisons of predicted and observed peculiar ve- 
locities. Previously, the tightest constraints on Pi were from 
the merged spiral and elliptical pecu liar velocity samples 
such as Mark III JWillick et al.lll997l) and SECat JZaroubil 
|200(J) as well as the SBF sample of lTonrv et all dl997l) . This 
work adds a result from local SNIa, a fourth independent 
data source of comparable statistical power. Recent compar- 
isons of predicted and observed peculiar velocities ('velocity- 
velocity'), including the result presented here, all yield re- 
sults consistent with a value of Pi — 0.5. 

Some of the earliest estimates of P were obtained by 
matching the gravity at the LG to the measured CMB 
dipole. While the LG has the most accurate observed CMB- 
frame velocity, a weakness of this method is that one needs 



to integrate the density field over all space to obtain the pre- 
dicted gravity at the LG. This contrasts with the velocity- 
velocity comparison performed above in which large-scale 
contributions to the predicted peculiar velocities either drop 
out of the analysis (if the fits are performed in the LG frame) 
or can be fitted independently of P (if the fits are performed 
in the CMB-frame) . This degeneracy cannot be broken when 
using the LG alone as one would be attempting to fit 4 pa- 
rameters (P and three components of an external dipole) 
to 3 degrees of freedom (the Cartesian components of the 
LG's CMB-frame motion). Consequently, in order to apply 
this method one needs either a deep, full-sky redshift survey 
(so that the external dipole is known to be zero) or, failing 
that, accurate estimates of the uncertainties arising from 
shot noise at large distances and from incompleteness in the 
ZoA. As an example of the latter. Irludson et alJ i2004h have 
sugges ted, based on the "Behi nd the Plane" extension of the 
PSCz iSaunders et al . 2000a), that additional structure in 
the ZoA beyond 100 /i -1 Mpc may increase the PSCz dipole 
by ~ (170±85) km s _1 . Until these issues are fully resolved, 
P determinations by this method remain subject to larger 
systematic errors than velocity- velocity comparisons. 

It is also possible to estimate P by comparing the den- 
sity field inferred from redshift surveys, to the mass density 
field constructed from peculiar velocity data. Such compar- 
isons are difficult because the mass density field is based on 
an inversion of sparse and noisy peculiar velocity samples. 
The POTENT rec onstruction of the Mark III catalogue by 
ISiead et al.1 lll998fr yields Pi = 0.89 ± 0.12. However, a new 
improved inversi on method based on an unbiased variant of 
the Weiner filter dZaroubi et all2002h finds Pi = 0.57±0.12, 
in good agreement with the velocity-velocity results. 

An alternative estimate of Pi can be obtained from 
other independent analyses not directly based on peculiar 
motion studies. One noteworthy route is via the combina- 
tion of parameters: Q^'as, where ag is the rms amplitude 
of mass fluctuations, 8 m , averaged within a top-hat sphere 
of 8 /i -1 Mpc radius. This combination may be related to 
Pi by the dependence of as,i, the number density fluctua- 
tion of IRAS galaxies, on the bias parameter bi. Since we 
are assuming linear biasing, Si = bi8 m and it follows that 
&s,i = bias. We can thus write: 

Pi = —r— = 

oi as,i 

ISpergel et, al.1 l)2003fi have used data from WMAP and other 
CMB and non-CMB sources to derive a value of Q^ag = 
0.38ln'n^. By directly integr ating the PSCz power spectrum 
iHamilton fc Tegmarld l|2002|) found a 8 .i = 0.80 ±0.05. Com- 
bining these two results gives Pi — 0.48 ± 0.06. The good 
agreement of the results from all these methods suggests 
that Pi is now known at the 10% level. 
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Table 2. (3j dependency on sampling 



Sample 




No. SNIa 


ft 




h' 1 Mpc < distance < 150 h~ 


1 Mpc 


98 


0.55 ± 0.06 


98 


h' 1 Mpc < distance < 30 h' 1 


Mpc 


31 


0.55 ± 0.07 


26 


30 h' 1 Mpc < distance < 150 h 


_1 Mpc 


67 


0.54 ±0.10 


74 


20 h- 1 Mpc < distance < 150 h 


_1 Mpc 


80 


0.55 ± 0.07 


84 


40 h' 1 Mpc < distance < 150 h 


_1 Mpc 


60 


0.49 ± 0.13 


67 


h' 1 Mpc < distance < 100 h~ 


1 Mpc 


85 


0.58 ± 0.06 


78 


h-' 1 Mpc < distance < 125 h~ 


1 Mpc 


90 


0.56 ± 0.06 


84 


No A v cull 




107 


0.50 ±0.08 


141 


A v < 0.5 




80 


0.57 ±0.06 


79 


A v < 0.3 




58 


0.57 ±0.08 


57 


CMB frame + dipolc 




98 


0.48 ± 0.09 


98 



Table 3. Recent determinations of j3i from velocity- velocity comparisons 



Comparison 



ft 



Reference 



Mark III vs. IRAS 1.2 Jy 
SNIa vs. IRAS 1.2 Jy 
SBF us. IRAS 1.2 Jy 
Mark III vs. IRAS 1.2 Jy 
Mark III vs. PSCz 
ENEAR vs. PSCz 
SFI vs. PSCz 
SEcat vs. PSCz 
SNIa vs. PSCz 



0.50 ±0.10 
0.40 ±0.15 
42+ ' 10 
0.50 ± 0.04 
0.60 ±0.10 
0.50 ±0.10 
0.42 ± 0.04 
0.51 ±0.06 
0.55 ± 0.06 



Davis et al. ( 1996*1 
Riess et al. (1997) 
Blakeslee et al. (19991 
Willick fc Strauss (19981 
Saunders et al. (1999 1 
Nusser et al. (2001) 
Branchini et al. (20011 
Zaroubi et al. (20021 
This Study 
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